Tunneling density of states of high Tc superconductors 
d-wave BCS model vs. SU{2) slave boson model 
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Motivated by recent experimental measurements of the tunneling characteristics of high Tc mate- 
rials using scanning tunneling spectroscopy, we have calculated the IV and differential conductance 
curves in the superconducting state at zero temperature. Comparing the two results obtained via 
BCS-like d-wave pairing and the SU{2) slave boson approach, we find that the slave-boson model 
can explain the asymmetric background observed in experiments. The slave-boson model also pre- 
dicts that the height of the conductance peak relative to the background is proportional to the hole 
doping concentration x, at least for under-doped samples. We also observe the absence of the van 
Hove singularity, and comment on possible implications. 
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PACS numbers; 74.20.Mn, 74.25.Jb, 74.50. -hr 
I. INTRODUCTION 

Tunneling spectroscopy has been one of the funda- 
mental tools in studying the superconducting state of 
the high Tc materials. In recent years it has been pos- 
sible to use the scanning tunneling microscope (STM) 
to perform reproducible experiments on single crystal 
Bi2Sr2CaCu208+5 cleaved in ultra high vacuum ||,|^. In 
contrast to photoeniission experiments, which are local 
probes in wave-vector space, the STM is local in real 
space. Thus it does not provide any information that 
depends on momentum, however, it has much higher en- 
ergy resolution. The density of states (DOS) obtained 
from the ^ curve is a direct fingerprint of the single 
particle microscopic physics in the material and hence of 
some importance in constraining the possible theoretical 
models explaining the elusive high Tc physics. 

The STM spectra for the superconducting phase ex- 
hibit unusual structure in the DOS when viewed in the 
light of BCS even if effects such as energy dependence 
of the normal state density of states of sample and/or 
tip, existence of bandwidth cutoffs, unequal work func- 
tions of tip and sample and energy-dependent transmis- 
sion probabilities are included O] . One notable feature in 
the tunneling spectra is the asymmetric background with 
an enhancement for hole-tunneling into the sample. This 
feature is reproduced very well in our calculation based 
on the slave-boson theory. The slave-boson model also 
predicts that the strength of the background in the tun- 
neling spectra does not scale with the doping x while the 
sharp conductance peak scales linearly with x. Thus by 
measuring the relative strength of the background and 
the sharp conductance peak as a function of doping x, 
one can distinguish the BCS theory and the slave boson 
theory experimentally. 

We also observe the absence of the van Hove singular- 
ity, and comment on possible implications. One possible 



implication is particularly intriguing and consistent with 
Photoemission results. That is the quasiparticles have a 
long lifetime, r > 0.5/meV, below the superconducting 
gap, and a very short lifetime, r < 0.05/meV, (spin- 
charge separation) above the superconducting gap. 

II. D-WAVE BCS 



The differential conductance ^ displays in the sim- 
plest case of constant DOS in the tip and energy inde- 
pendent transition probability the single electron DOS 
in the sample. This reflects the ability of the material to 
accommodate an extra electron or hole depending on the 
sample bias. Within the tunneling Hamiltonian formal- 
ism the tunneling current is given by B 

JT = iireT^ ^ / c?^ i^L- {OJ, p)Ar+ {uj + V,k) 
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where the Aj^s and A^s are the spectral functions for 
the single electron Green's functions in the tip (L) and 
sample (R) respectively. V denotes the bias of the sam- 
ple with respect to the tip and F is the tunneling matrix 
element assumed independent of energy. Notice that pos- 
itive V corresponds to e~ tunneling into the sample. 

For a free fermion system which we suppose to rep- 
resent the tip material we have the standard form for 
the spectral function at zero temperature Alj^{ll!,p) — 
Q{uj)S{lj — ^p) and Al-{u>,p) = Q{—Lu)d{uj — ^p), were 
S.p — £p — Ep denotes the particle spectrum in the tip 
with Ep the Fermi energy. Q{lo) is the Heaviside step 
function where to is measured with respect to the Fermi 
energy. 

For the sample, we first consider the following single 
particle spectral distribution at zero temperature. 



AR+iuj,p) = e{uj)u'^{p)S{uj - Ep) 
Ar^[lo,p) ^ Q{-w)v\p)6{uj + Ep) 



(2) 



where u^{p) = Ul + ^) and ^^(p) == Ul - ^) are the 



BCS coherence factors Ep ~ J^"^ + A^ is the dispersion 
relation of the quasiparticles in the superconducting state 



and S,p 



Ep with ep the dispersion relation in the 



normal state. Here Ap denotes the gap function which 
is taken to have a d-wave symmetry in reciprocal space. 
Within the above approximation to the spectral functions 
we obtain the following expressions for the single particle 
tunneling current 

3T\v>o^^enT^N{EF) Y. ^'(^) 

k,E{k)<V 

Jt\v<o = -4enr^N{EF) ^ ^^(fc) 

k,E{k)<\V\ 

In Fig. (||) we have plotted the resulting differential 
conductance curve. With 



£.k — to + ti(cos kx + cos ky) + t2(cos 
+ i-? fcos 2k^ + cos 2fc„ + Ia fee 



is (cos 2kx + cos 2ky + t^ (cos 2kx cos ky 
+ cos 2ky cos kx) + ^5 cos 2kx cos 2ky — Ep 

as the dispersion relation for the quasiparticles in the 
ab - plane of the sample and the matrix elements cho- 
sen as follows [to,...,h] = [0.1305,-0.2976,0.1636,-0.026,- 
0.0559,0.051] in eV, this is a tight binding fit to angle 
resolved photo-emission spectroscopy (ARPES) measure- 
ments performed by Norman ei ai B. Ep has been ad- 
justed to yield 10% hole doping and A^ — Ao(cosfc:E — 
cos ky) { with Aq = 22meV) has the d-wave k-space sym- 
metry mentioned previously. 

With the hole doping at 10%, the van Hove singularity, 
present in the band structure, ends up on the hole side 
of the DOS very close to the Fermi energy. The fact that 
the van Hove singularity is close to the Fermi surface and 
hence should show up in the low energy single particle 
physics can be seen nicely in Fig. (^ in the form of the 
double peak structure. The coherence factors which mix 
particle and hole density of states lead to the van Hove 
singularity also showing up on the particle side of the ^ 
curve albeit with much smaller amplitude. 



III. SLAVE BOSONS 

Next we consider the tunneling problem in the light 
of the SU{2) slave boson theory of Wen and Lee [p|. It 
is commonly believed that the simplest model that in- 
corporates the strong correlation physics relevant for the 
high Tc cuprates is the t-J model. Due to the strong 
on-site Coulomb repulsion energy the doubly occupied 
states should not contribute to the low energy effective 
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FIG. 1. The two peaks symmetrically located in height 

and energy with respect to zero bias are the usual peaks aris- 
ing from the gap structure in the superconducting density of 
states. The two outer peaks are the remnants of the van Hove 
singularity. The wiggles in the background are due to the dis- 
creteness of k-space when performing numerical calculation. 
The discreteness gets amplified by the derivative taken in ob- 
taining ^ from the tunneling current (pi). The resolution in 
voltage is 3meV 



theory. Within the SU{2) approach this constraint is im- 
plemented via the introduction of a slave boson doublet. 
The physical electron operator can then be written as an 
SU{2) singlet. Within this representation, the mean-field 
electron propagator is given by the product of the boson 
and the fermion propagators and was calculated in M . 

For our purpose we only need the T = spectral func- 
tions which can be read off from the expression for the 
Green's function as 

An+{u;, k) = Q{u:)[^u)(k)5{uj - Ef {k))] 
for the particle part 

Ar_{uj, k) = Q{-uj)V^v){k)5{Lo + Ef{k)) 

+ 27V l^t"'''^'^ ^ ^)"/(<?) + «b('? - k)vf{q)Y 
1 

X S{uj + Ef {q) + Et{q - k)) (3) 

+ 2iV 5I["&<^* ^ k)vf{q) - Vb{q - k)uf{q)f 



X S{LU + Ef{q)+El{q-k))^ 



for the hole part of the spectrum. Here N de- 
notes the number of sites and x is the hole dop- 
ing concentration. The remaining variables are de- 
fined as follows: M/,b(fc) = ^^y 1 + fj/.tp «/,b(fc) = 

and £:^(fc) = ±\/(e''(fc))2 -I- (A''(fc))2 - ^f,, where 
e-^{k),A^{k) are the fermion dispersion and gap func- 
tion respectively e^{k),A^{k) are the dispersions of the 
bosons and /x;, is the boson chemical potential. 



Notice that besides a coherent part for the spectral 
functions which resembles the form of the BCS spectral 
weight (g) albeit scaled by a factor of |- there is also 
an added incoherent contribution to the hole part of the 
spectral function (||). 

To calculate the tunneling current we have used the fit 
to ARPES measurements as dispersion for the fermions 
(spinous) and a nearest neighbor tight binding dispersion 

e^ = —2t^{coskx + cos ky) 

for the bosonic degrees of freedom (holons) with t^ the 
hopping matrix element for the holons. It is important 
that we match the fermionic band structure with the 
ARPES measurements since the fermions have a bigger 
band mass and hence determine the dispersion relation 
seen in ARPES |^,^ The electrons measured in those ex- 
periments are thought of (within spin-charge separating 
models) as bound states of the heavy spin degrees of free- 
dom and the light charge degrees of freedom. Since we are 
interested in the low energy effective theory the details 
of the broad dispersion for the holons (charge degrees of 
freedom) are not crucial and hence have been chosen as 
simple as possible. Furthermore to arrive at equations 
(^ we assumed boson condensation of the holons. 

With the above expressions for the spectral functions 
in the sample we can calculate the tunneling current us- 
ing equation dl^) as 

X 
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JT\v<o = -4e7rT^NiEF) ^ ^v^fik) (4) 

k,E{k)<\V\ 



2N 
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k,q 



where Ek, = Ek^^Hq) + E^q - k)) - e{Ef{q) + 
Ei{q-k)-\V\)]. 

When comparing the two ^ curves Fig. (|l|) and (||) 
one can see how the lowest energy physics is virtually 
identical. However on energy scales bigger than 4Ao 
a marked asymmetry in the background of the SU{2) 
model 4^ shows up with an increase in the hole tunnel- 
ing spectral weight. 

The inset of Fig. depicts the incoherent contribution 
to the hole tunneling spectrum separately. Scaling the 
height of Fig. by | and adding the incoherent hole 
contribution results in Fig. 0. The increase in the hole 
tunneling spectral weight arises due to the fact that re- 
moving an electron from the sample requires the recom- 
bination of the spin and charge degrees of freedom into 
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FIG. 2. 



Sample Bias (mV) 

re have set t{/2t'' 



Here we have set t(/2i'' = 1/2. 
A-' = Aq (cos kx — cos ky) and A' = Ao(cos kx — cos ky) with 
A^/Af) = 1/2 where t{ = - 297.6 meV from the ARPES 
fit and Aq — 22 meV. Notice that above ratios correspond 
to J/t — 1/2 within the t-J model. The solid curve is our 
slave-boson result. The wiggles in the background are more 
pronounced here as compared with Fig. (fll) due to smaller 
resolution in k-space in calculating the convolution integrals 
(y). Here the resolution in 1/ is 5 meV. The dots are the 
experimental result from Ref. [2]. The inset shows the inco- 
herent part of the hole tunneling spectral weight. Notice that 
the exact shape of this curve should not be taken too literally 
- see discussion on boson band structure in the main text. 



a single entity. This yields a mixing in of the higher en- 
ergy holon dispersion whose detailed form is not known 
within the effective low energy theory. 

Another feature so far not discussed is the scaling with 
the hole doping of the conductance peak corresponding 
to electrons tunneling into the sample. Comparing equa- 
tions (^) and (|j) we see that within the SU{2) model, 
the peak height scales linearly with x whereas there is 
no dependence of the peak height on doping within the 
d-wave approach. The doping dependence within SU{2) 
arises from the reduction of the overlap of the electron in 
the tip with the quasiparticle (as a bound state of holons 
and spinous) in the sample which crudely speaking means 
that an electron can only enter the sample on empty sites 
and then 'decay' into its constituent parts. The only dop- 
ing dependence within the d-wave approach arises due to 
the chemical potential which dictates the separation of 
the double peak structure but not its height. The lin- 
ear scaling with x within the SU{2) slave boson mean- 
field theory discussed here should be taken more as a 



qualitative than exact quantitative prediction, since it is 
a mean-field result. Recent photoemission experiments 
|g[ by Ding et al observed a linear x dependence of the 
quasiparticle peak, which fits the mean-field result of the 
SU{2) theory very well. 

Thus, at the mean-field level we have found qualitative 
different behaviors with regards to the x dependence of 
the hole tunneling background and the electron tunnel- 
ing peak within d-wave BCS and the SU{2) slave boson 
theory. It is this difference in doping dependence, which 
should be experimentally testable and hence yield to a 
feature distinguishing between the two models. 

Furthermore notice that the DOS contains singularities 
(for both BCS and the slave boson model) at the electron 
tunneling peak. The curvature of the measured dl/dV 
curve at these peaks should give us an upper bound on 
the quasiparticle decay rate at the energy scale of the 
superconducting gap. Based on new experimental data 
by Pan et al 0, the quasiparticle decay rate can be as 
small as a few meV even for quasiparticles with energy 
as high as AQmeV . This is very different from the normal 
state where the quasiparticle decay rate is of the same or- 
der of magnitude as the quasiparticle energy. We would 
also like to remark that according to the photoemission 
and tunneling results for underdoped samples, the quasi- 
particle peak (with a width of order T) disappears com- 
pletely above Tc while the gap remains at (0,7r). Based 
on the slave-boson theory, the sharp electron-tunneling 
peak arises due to the condensation of the holons (whose 
weight is proportional to a; at T = 0). As T approaches 
Tc^ the fraction of the condensed holons vanishes. If we 
assume that the holons are very incoherent above T^ 
we can conclude that the sharp electron tunneling peak 
should disappear above Tc- This picture from the slave- 
boson model is completely consistent with the observed 
results from photoemission experiments. 

Another point to make here is about the van Hove 
singularity. Samples with small superconducting gaps 
(A ~ 25meV) show a double peak structure in the 
tunneling dl/dV curve, and the double peak structure 
crosses over into a single peak for large superconduct- 
ing gaps. At first sight, one might guess that the double 
peak structure is due to the van Hove singularity. How- 
ever, after comparing the experimental lineshape with 
the theoretical lineshape, we conclude that the double 
peak cannot arise due to the van Hove singularity. This 
is because the experimental peaks at higher bias are quite 
symmetric, while the peaks from the van Hove singularity 
are very asymmetric (the peak at the hole side is much 
stronger than the peak on the electron side). Based on 
the dispersion obtained from the fitting of the ARPES 
measurements, the van Hove singularity should show up 
even for samples with larger gaps (A — 50meV). How- 
ever, experimentally, one fails to see the van Hove singu- 
larity even when the gap is as small as 20meV. H This 
seems to suggest that quasiparticles have very short life- 



time (spin-charge separation) above the superconducting 
gap, and hence the van Hove singularity cannot be ob- 
served. This leaves us with the question of where the 
double-peak feature comes from if there are no well de- 
fined quasiparticles above the superconducting gap? We 
hope that the double-peak structure may give us some 
hints on how coherent quasiparticles emerge in the su- 
perconducting state from the incoherent normal state. 

Finally we would like to point out that the results 
of this paper are obtained from a mean-field calculation 
within the slave-boson theory. One naturally questions 
the reliability of the mean-field result and how much of 
our result remains valid after the gauge and other fluc- 
tuations are included. One of our main results is the 
explanation of the asymmetric tunneling background. It 
can be traced back to the strong on-site repulsion. It is 
much easier to remove an electron than to add an elec- 
tron and create a doubly occupied site. We believe this 
result is robust and will survive the fluctuations around 
the mean-field state. The second main result is that the 
weight of the coherent quasiparticle tunneling peak is 
proportional to the doping x. After including the fiuctu- 
ations, we believe that the weight of the coherent peak 
should have a similar doping dependence. This is because 
the coherent peak comes from the quasiparticle which is 
a bound state of a spinon and a holon. However, the de- 
tailed dependence may be of a more general form a;^"*"" 
{ie the fluctuations may correct the exponent). pOJ 

We would like to acknowledge helpful discussions with 
J.C.Davis, D.A Ivanov, P.A. Lee and S.H. Pan. This 
work is supported by NSF Grant No. DMR-97-14198 
and by NSF-MRSEC Grant No. DMR-98-08941. 



[1] C. Renner et al, Phys. Rev. Lett. 80, 149 (1998). 

[2] S.H. Pan et al, Nature 403, 746 (2000) 

[3] J.E. Hirsch (1999), preprint (cond-mat/9901059). 

[4] G.D. Mahan, Many-Particle Physics (Plenum Press, New 
York, 1990) 

[5] M.R. Norman et al. Nature 392, 1571 (1998). 

[6] X.G. Wen and P.A. Lee, Phys. Rev. Lett. 76, 503 (1996). 

[7] P.A. Lee, N. Nagaosa, T.K. Ng and X.-G. Wen, 
Phys. Rev. B 57, 6003 (1998). 

[8] H. Ding et al, |3ond-mat/0006143 . 

[9] R.B. Laughlin, Proceedings of the Inauguration Confer- 
ence of the Asia-Pacific Center for Theoretical Physics, 
Seoul National University, ed. by Y.M. Cho, J.B. Hong, 
and C.N. Yang(World Sci., Singapore, 1998) 
[10] Walter Rantner and Xiao-Gang Wen, to appear. 



